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We consider the Sinai model describing a particle diffusing in a one dimensional random force 
field. As shown by Golosov, this model exhibits a strong localization phenomenon for the thermal 
packet: all thermal trajectories starting from the same initial condition in the same sample remain 
within a finite distance of each other even in the limit of infinite time. More precisely, he has proved 
that the disorder average Pt(y) of the distribution of the relative distance y = x — m{i) with respect 
to the (disorder-dependent) most probable position m(t), converges in the limit t — > oo, towards 
a distribution Pg(u) defined as a functional of two independent Bessel processes. In this paper, 
we revisit this question of the localization of the thermal packet. We first generalize the result of 
Golosov by computing explicitly the joint distribution Poo(y, u) of relative position y = x(t) — m(t) 
and relative energy u = U(x(t)) — U(m(t)) for the thermal packet. Next, we compute the localizations 
parameters Yk, representing the disorder-averaged probabilities that k particles of the thermal packet 
are at the same place in the infinite-time limit, and the correlation function C(l) representing the 
disorder-averaged probability density that two particles of the thermal packet are at a distance 
/ from each other. We moreover prove that our results for Yk and C(l) exactly coincide with 
the thermodynamic limit L — > oo of the analog quantities computed for independent particles at 
equilibrium in a finite sample of length L. So even if the Sinai dynamics on the infinite line is always 
out-of-equilibrium since it consists in jumps in deeper and deeper wells, the particles of the same 
thermal packet can nevertheless be considered asymptotically as if they were at thermal equilibrium 
in a Brownian potential. Finally, we discuss the properties of the finite-time metastable states that 
are responsible for the localization phenomenon and compare with the general theory of metastable 
states in glassy systems, in particular as a test of the Edwards conjecture. 



I. INTRODUCTION 

The Sinai model jij of a particle diffusing in a one di- 
mensional quenched random force field is one of the sim- 
plest example of a model with quenched randomness. Its 
continuum version is defined by the Langevin equation 



dx(t) 
dt 



-U'[x(t)]+rj(t) 



(1) 



where ry(t) is the thermal noise, with correlation < 
r](t)r)(t') >= 2T5(t — t'), and where the random potential 
U(x) is a Brownian motion presenting the correlations : 



(U{x) - U{x')) 2 = 2a\x-x'\ 



(2) 



As a result, the Sinai diffusion exhibits a non trivial ultra- 
slow logarithmic behavior, the walker typically moving as 
x ~ (lnt) 2 - Although this model has been much stud- 
ied QH^, the known analytical results mainly concern 
the rescaled variable X — t^St? * , and its distribution 
over the disorder realizations, known as the Kesten dis- 
tribution. However another important issue concerns the 
thermal distribution of the position in a given sample. 

Golosov g] has discovered the important phenomenon 
of localization in the sense that all thermal trajectories 
starting from the same initial condition remain within 



a finite distance of each other even in the limit of infi- 
nite time! More precisely, he has proved that there ex- 
ists a process m(t), independent of the thermal noise 77, 
such that the distribution Pt(y) of the relative distance 
y = x — m (t) averaged over the realizations of the disor- 
der converges in the limit t — > 00 towards a probability 
distribution Pg{v) defined as the following functional 



Pa(y) 



-r(|v|) 



f °° dte-rW + J Q °° dte- 



p(t) 



(3) 



{r,p} 



where << .. >> denotes the average over the two in- 
dependent Bessel processes (i.e. the radial parts of free 
Brownian motions in three dimension) r(t) and p{t) start- 
ing at r[t = 0) = = p(t = 0). (The explicit compu- 
tation of this functional is done in Appendix ^ of the 
present paper.) 

However, the existence of this limit distribution 
doesn't imply that the moments of the random vari- 
able y remain finite in the limit of infinite time. And 
indeed, all the integer moments of the relative dis- 
tance (x— < x(t) >) to the thermally averaged position 
< x(t) > diverge in the infinite-time limit, with the fol- 
lowing leading divergence |t],^| 



T 



<(x-<x(t) >T>~—{T\ntY 



(4) 
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where < .. > denotes the thermal average over {r)(t)} 
and where — denotes the disorder average over the ran- 
dom potential {U(x)}. This happens because the decay 
of the distribution of z = (x— < x(t) >) is algebraic at 
large distance as l/|z| 3 / 2 For n = 2, the behavior 

(|J) has been measured numerically in || . 

Other quantities characterizing the localization of the 
thermal packet are the localization parameters Yk (t) rep- 
resenting the disorder-averages of the probabilities that 
k independent particles in the same sample starting from 
the same initial condition are at the same place at time t. 
In a given environment U (x) and for a given initial con- 
dition xq, the probability distribution over the thermal 
noise 



P(x,t\x ,0) =< S (x - X{ VtUiX0 y(t)) > 



(5) 



(where X{ v jj xo y(t) is the solution of the Langevin equa- 
tion (0)) satisfies the Fokker-Planck equation 



d t P(x, t\xo,0)= -H FP P(x, t\x ,0)) 
H FP =-d x (Td x + U'(x)) 

and the initial condition P(x,t — > 0|xo,0) — 
So the localization parameters read 



Y k (t) 



+00 



dx[P(x,t\x o ,0)} 



(6) 
(7) 

S(x - x Q ). 



(8) 



In ||, the parameter Yzit) has been measured numeri- 
cally for a version of the discrete Sinai model in a semi- 
infinite geometry with binary distribution of the random 
forces. This simulation shows that Yi{t) converges at 
long time towards a finite value 12(00), which decays as 
T increases (since the temperature broadens the distri- 
bution of the thermal packet). 

A generalization of l^(^) is the correlation function 
C(l, t) representing the disorder-average of the probabil- 
ity that two independent particles in the same sample 
starting from the same initial condition are at a distance 
I from each other at time t 



C(l,t) 



dx[P(x,t\x o ,0)P(x + l,t\x o ,0)] (9) 



which is normalized to J Q dlC%{l,t) = 1. 

Another question related with the distribution of the 
thermal packet is the dynamics of a given particle be- 
tween two times t w and (t w + r) in the 'quasi-equilibrium 
regime' t w — > 00 with finite r I. It was conjectured 
and checked numerically in Q that the disorder aver- 
aged probability Q(z,t) = limt w _ >00 Q(z, t w + r,t w ) of 
the relative displacement z = x{t) — x(t w ) for the Sinai 
model on the infinite line was the same as Q e q( z , T ) ob- 
tained as the thermodynamic limit L — > 00 of Ql{z,t) 
characterizing the equilibrium dynamics in a finite sam- 
ple of length L. In particular, in the large r limit, one 
should recover the statics with two independent parti- 
cles at Boltzmann equilibrium M. If these assumptions 



are true, this means that for the particles of the thermal 
packets, we should also have a correspondence with the 
Boltzmann distribution in a Brownian potential on finite 
sample of length L 



e -HU{x) 

p e r q (x) — —f 

tfdye-mv) 



(10) 



More precisely, it is interesting to compare Yfc(oo) (j|) 
and C(Z,oo) (|^) with the thermodynamic limit of their 
statics counterparts 



Y? = lim / dx[p e ^ (x)] k 



(11) 



L — ►oo 



C^{1)= Urn / dx I dyp eq (x)pl q (y) (12) 



Some statistical properties of the Boltzmann distribution 
( |l0| ) have already been studied in |pl]-p^[| . 

In this paper, we reconsider this question of the local- 
ization of the thermal packet from the point of view of 
the real-space RG analysis detailed in Q|. Within this 
renormalization picture, at time t, any particle starting 
from an initial condition belonging to a renormalized val- 
ley will be typically at time t around the minimum m(t) 
of this renormalized valley. To study the distribution of 
a thermal packet, a first step is to consider that the par- 
ticles of the packet are at Boltzmann equilibrium within 
the renormalized valley they belong to at time t. This is 
only an approximation at finite time, since there are also 
additional out-of equilibrium situations for the thermal 
packet 0. However in the limit of infinite time, these 
out-of equilibrium situations have vanishing probability 
, and the joint distribution Poo (y, u) of relative position 
y = x(t)—m(t) and relative energy u = U(x(t)) — U(m(t)) 
corresponds to an average of Boltzmann distribution over 
infinitely deep Brownian valleys. We will compute ex- 
plicitly P 00 (y,u) by a path-integral method. We use the 
same method to compute the Yfc(oo) parameters (H) and 
the correlation function C(l) (H ). On the other hand, 
we compute Y^ q and C eq (l) ( |12"| ) and find that they in- 
deed coincide with Yfc(oo) and C(/,oo). This shows that 
the ensemble of infinitely-deep valleys gives the same re- 
sults for the quantities mentioned above as the thermo- 
dynamic limit of the ensemble of finite-size valleys, so 
that quasi-equilibrium in Sinai diffusion and equilibrium 
in a Brownian potential are equivalent. 

This approach to the localization phenomenon allows 
us to study in details the disorder-dependent structure 
of low-energy eigenstates of the Fokker-Planck operator. 
Our results are consistent with the features discussed in 
the Appendix B of |l4| for the related model of one- 
dimensional random-hopping hamiltonian for fermions. 

Finally, it is instructive to recast Sinai diffusion into 
the general theory of glassy systems. Indeed, in the stud- 
ies on slowly relaxing systems such as glasses, granular 
media or disordered spin models, it is natural to sepa- 
rate the dynamics in two parts : there are 'fast' degrees 
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of freedom which rapidly reach local quasi-equilibrium 
plus a slow non-equilibrium part. At a given long-time 
t, the fast motion covers a region of phase-space which 
can be defined as a metastable state associated to time 
t [p^[ . In Sinai diffusion, this picture directly applies : 
the metastable states are the renormalized valleys, within 
which there is a local Boltzmann equilibrium. Moreover, 
we obtain that the metastable states satisfy all the prop- 
erties of the construction |l6| as summarized in |H| . In 
glassy systems, the Edwards ergodicity conjecture ]17f| 
consisting in computing dynamic quantities by taking flat 
averages over metastable states has given rise to a lot of 
recent studies p^^!q -pi|. Since Edwards conjecture is 
usually based on the assumption that all the basins of 
attraction of the various metastable states have the same 
size , it is of course a very strong hypothesis that can- 
not be true in general but only for special systems with 
special dynamics ]l8| , pl| . In Sinai diffusion with uniform 
initial condition, the size of the basin of attraction of a 
metastable state is given by the spatial length of a renor- 
malized valley : it is thus a random variable whose distri- 
bution is exactly known. As a consequence, Edwards con- 
jecture cannot be true in general. Nevertheless, within 
the RSRG approach , all one-time quantities are effec- 
tively computed by averages over all metastable states, 
but with a measure which is not flat, but depends on the 
quantities and on the properties of the basins of attrac- 
tion. So the RSRG approach of the Sinai model repre- 
sents the simplest example where the dynamical study of 
a glassy system can be faithfully replaced by an average 
over a set of well specified metastable states, with a well 
defined measure. However, for special quantities, Ed- 
wards conjecture can be recovered. For instance, in this 
paper, we compute the probability distribution Poo{y,u) 
of the thermal packet, the localization parameters Yfc(oo) 
(||) and the correlation function C(Z, oo) (^|) as flat aver- 
ages over infinitely deep wells. This is because in the 
infinitely deep valleys, the statistics of the lower-part of 
the Brownian valley is the same for all metastable states. 

The paper is organized as follows. In section 0, we 
explain within the Real-space Renormalization Group 
(RSRG) picture why the distribution of the thermal 
packet is asymptotically given by an average of Boltz- 
mann distribution over infinitely deep Brownian valleys. 
In section [II, we use a probabilistic path-integral method 
to compute explicit expressions for the joint probability 
distribution Poo(y, u). We use the same method to com- 
pute the probability distribution of the partition func- 
tion of an infinitely-deep valley (section [jv|) , the local- 
izations parameters Yfc(oo) (section [v]), and th e co rre- 
lation function C(Z,oo) (section |v|) . In section VII, we 
consider equilibrium functions in finite samples and com- 
pute Y^ q and C eq (l) ( |ll| ) which ar e foun d to coincide 
with Yfc(oo) and C(/,oo). In section VIII, we derive ex- 



appendices [A|, |b| and |c| contain technical details used in 
the text, whereas the appendix |d| shows that our result 
for Poo(y,u) after integration over u coincides with the 
explicit computation of the Golosov functional (^). 



II. REAL-SPACE RENORMALIZATION GROUP 
FOR THE SINAI DIFFUSION 



In this section, we briefly recall the principles of the 
Real-Space Renormalization Group approach to Sinai dif- 
fusion 0| with special emphasis on the successive levels 
of approximations. 



A. Effective dynamics at large times 



Recently we have proposed an approach, based on a 
real space renormalization group (RSRG) method, which 
allows us to obtain many exact results for the Sinai model 
HQ . The way to implement the RSRG is very direct: one 
decimates iteratively the smallest energy barrier in the 
system stopping when the time to surmount the small- 
est remaining barrier is of order the time scale of inter- 
est. Despite its approximate character, the RSRG yields 
for many quantities asymptotically exact results, because 
the iterated distribution of barriers grows infinitely wide. 
Indeed, the distribution of the rescaled barrier height 
rj = f p F converges towards the fixed point 



where 



P*(rj) = %) e -" 



T = Tint 



(13) 



(14) 



plicit expressions for the eigenfunctions of the Fokker- 
Planck operator. In section IX, we discuss the properties 
of metastable stables and compare with the general the- 
ory of metastable states in glassy systems. Finally, the 



is the renormalization scale associated to the time t. 

Within this renormalization picture, at time t, any 
particle starting from an initial condition belonging to 
a renormalized valley (F 1 ,F 2 ) will be typically at time 
t around the minimum m(t) of the valley. This sim- 
ple approximation, called "effective dynamics" in [Q, is 
sufficient to obtain exact expressions for many quanti- 
ties, such as for instance the Kesten distribution of the 
rescaled variable X — J^S*}. . 

However for other quantities, we have already obtained 
in j7j] that there are differences between the effective dy- 
namics and the real dynamics. For instance, persistence 
properties of the thermal average < x(t) > are well de- 
scribed by persistence properties of the effective dynam- 
ics consisting in jumping between valley bottoms but are 
very different from the persistence properties of a single 
walker @. 
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B. Boltzmann equilibrium within renormalized 
valleys 

To study the distribution of a thermal packet, we 
clearly need to go beyond the effective dynamics. A first 
step is to consider that the particles of the packet are at 
equilibrium within the renormalized valley they belong 
to at time t. More explicitly, this approximation which 
assumes that the walkers are at Gibbs equilibrium sep- 
arately in each renormalized valleys at scale T, can be 
written 



P(xt\x 0) -J^J-e-WWOnWOvrixo) 



(15) 



where the sum is over all the renormalized valleys Vr that 
are present in the system at the renormalization scale 
r = Tint, and where 9y(x) is the characteristic function 
of the valley V, i.e Ov(x) = 1 if x belongs to the valley 
and Oy(x) — otherwise. Zy = f v dxe~^ utyX ^ represents 
the Boltzmann normalization over the valley V. 

This Boltzmann equilibrium is clearly an excellent ap- 
proximation within the lower part of the valley, i.e. for 
the points that are at a finite potential above the min- 
imum of the valley, which have had plenty of time to 
equilibrate. However, it breaks down further away in the 
higher part of the valley, for the points that are at a po- 
tential of order T above the minimum of the valley, since 
these points need a time of order e@ r ~ t to equilibrate. 
However, since these points have a weight of order e~~^ r 
in the formula (|l5|), they do not play any role for the 
observables computed in this paper in the limit r — > oo. 

More importantly, the approximation ( |l5| ) breaks 
down whenever out-of-equilibrium situations occur for 
the thermal packet as we now explain. 



C. Out-of-equilibrium situations for a thermal 
packet 

In our previous work we have already described 
rare events where deviations from the effective dynamics 
show up. The most important rare events are of order 
1 /r and are of three types pictured on Fig 7 of Ref. |Q . 
In the events of type (a), there are two nearly degener- 
ate minima at thermal equilibrium separated by a barrier 
To < r. These events (a) are thus well taken into account 
by the Boltzmann equilibrium in each renormalized val- 
ley described in previous section. The rare events (b) 
where two tops are nearly degenerate are on the contrary 
completely out-of equilibrium, since the thermal packet 
will be split in two valleys that are not at equilibrium 
with each other. Finally the events (c) where the val- 
ley is being decimated are also out-of-equilibrium events, 
since the two renormalized valleys at T cannot be consid- 
ered at thermal equilibrium. All other rare events are of 
higher order, for instance the probability that the initial 



point is near a top, which will also produce an out-of- 
equilibrium splitting of the thermal packet, is of order 

i/r 2 . 



D. Conclusion 

As a conclusion, the expression (|l5|) is an excellent 
approximation for the thermal packets that are not in 
out-of-cquilibrium situations, and within the lower part 
of the renormalized valleys, i.e. for the points that are at 
a finite potential above the minima of the valleys. This 
approximation breaks down for the higher parts of the 
valleys, i.e. for the points that are at a potential of order 
r above the minima of the valleys, and whenever the ther- 
mal packet happens to be in an out-of-equilibrium situ- 
ation like the events (b) and (c) described above which 
appear with probability 1/r. For a detailed study of sys- 
tematic corrections to this approximation, and further 
results, we refer the reader to D.S. Fisher p3fl . 

From now on, in this paper, we will restrict our atten- 
tion to the approximation ( |l5| ) which becomes asymptoti- 
cally exact in the limit T — + oo. Indeed, in this limit, only 
the points that are at a finite potential above the min- 
ima of the valleys keep a finite weight and all the out-of 
equilibrium situations have a vanishing probability in the 
limit r — > oo. As a consequence, the Boltzmann distri- 
bution over an infinite-deep valley exactly describes the 
asymptotic dispersion of a thermal packet. More explic- 
itly, the disorder average Poo(y, u) of the joint probability 
distribution of the relative position y = x(t) — m{t) and 
the relative energy u — U{x(t)) — U(m(t)) with respect 
to the minimum [m(t), U(m(t)] of the valley is given by 



e-0 u 6(u-Ui(\y\)) 



Jo°° dxe ~ 



/3£/i(a0 



+ Jo°° dxe ~ 



(16) 



{Ui,U 2 } 



where the average (...) is over two independent Brownian 
trajectories U\(x) and U2(x) forming an infinitely-deep 
well, i.e. a renormalized valley in the limit T — > oo. We 
note here that the minimum m(t) of the valley represents 
the most probable position in each sample (i.e. it is the 
point where the probability is the biggest), but not the 
thermally averaged position < x(t) >. 

Using, the same notations, we obtain the infinite-time 
limit of the localizations parameters ife(oo) (@) 

/■+0O 

Yfc(oo) = / dy 



and of the correlation function (m 



(17) 
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C(Z,oo)=4/ dy[P 00 (y)P 00 (y + l)} + 2 / dy^l - y)P 00 (y)} 
Jo Jo 
r°° I -pu^-pu^y+i) \ pi I 

= 4 / dy ( 5- ) + 2 / dy 

Jo \(f™dxe-Wi(x) + f™dxe-P u ^)) 2 / ' 



o 



(/ + °° dse-^O") + J +o ° dxe-W 2 (x)y 



(18) 



We now turn to the explicit computation of these ex- 
pressions. 



III. PROBABILITY DISTRIBUTION FOR THE 
THERMAL PACKET 

A. Expression of P oa {y,u) in terms of path-integrals 

To define more precisely the average (...) ( |l6| ) over 
two independent Brownian trajectories U\(x) and ^(x) 
forming an infinitely- deep well, let Vr be the set of Brow- 
nian paths {U(x > 0)} starting at U(0) = + in the 
presence of absorbing boundaries at and T, and that 
are conditioned to finish at U = T and not at U = 0. 
The formula ([iff) is thus defined as 



Poo{y,u) = lim 
r— »oo 



-/3m 



*(u-ffi(|l/|) 



,/ ' r dxe-^ x ) dxe-Wii*), 



(19) 



where U\{x) and U2{x) are two independent Brownian 
trajectories belonging to Vr and where l^ and l^ are 
the random times where U\{x) and U%{x) respectively 
first hit x = r where they are killed. Since the expres- 
sion is symmetric in y — > — y, we will assume y > from 
now on. 

To separate the averages over U\ and U2, it is conve- 
nient to exponentiate the denominator to get 



Poo 

(y,u)= lim R r (q)S r (y,u,q) 

T— >oo 



where 



R r {q) = le- q ^ 



(20) 



(21) 



{U2} 



and 



Sr{y,u,q) = 



(22) 



{Ui} 



F [0X] (uJ\u ) = / VU(x) 

JU(0)=u o 

e-^Ii^r-^U^-^e^iU} (23) 

where 0[o, r]{U(x)} means that there are absorbing 
boundaries at U = and at U = T. The explicit com- 
putation of this path-integral is done in Appendix ^ and 
yields the final result (B18). 

To compute the quantity (|2l]), we need to consider the 
path-integral (|2^) going from the initial potential uo = e 
to the final potential u = T — e in the limit e — > and 

to sum over the random time / representing the random 
(2) 

time ; where U2(x) first hit x — V where it is killed. 
So we have 



,+00 1 
R T (q)=Af(T) / dl]im^F [0tr] (T-e,l\e) 
Jo e 



(24) 



up to a normalization Af(T) that ensures R r(q = 0) = 1. 
The result for Rr(q) is given in equation ( B25 ) of Ap- 
pendix H which yields in the limit r — > 00 



Roo(q) 



1 



Hi 



(25) 



Similarly, to compute (B2), we need to compose two 
path-integrals of type (|2^), the first one going from the 
initial potential uq — e to the final potential u in a time 
y, and the second one going from the initial potential u 
to the final potential T — e in a time I representing the 
difference (Z r — y) that we have to sum over, so that we 
have 



S r (y,u,q)=M(T)e-^ 



+OO 



dl 



lim —F[o,r] ( r ~ e > l \ u ) F [o,r](u, y\e) 



(26) 



The Laplace transform wit h resp ect to y of this expres- 
sion is given in equation (B30) of Appendix |b] which 
yields in the limit T — > 00 



Soo(p,u,q) = 
I u (se~^ 



K (se 



dye py S r (y,u,q) 
v *o(«) 



-/3u 



Pa 



Io{s) 



I (se ^) 



(27) 



We now define the path-integral 



where 
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V cr 

1 /£ 

f3\ a 



(28) 
(29) 



B. Final result for the joint probability distribution 

Poo(y,u) 



probability of return to the origin of a random walk [pL[7j . 
Indeed, the Boltzmann distribution in a renormalized val- 
ley typically decays at large distance as e - ^^ 5 ^. How- 
ever, in rare configurations where the random potential 
U (y) happens to be near the origin U(0) at y, which hap- 
pens with probability l/j/ 3 / 2 , then the Boltzmann distri- 
bution has a weight of order 1 at y. These rare configu- 
rations entirely dominates the disorder-average for large 
y and are responsible of the power-law decay . 



Using the results ( |25|) and ([27J) , the Laplace transform 
of the probability distribution P^ {y,u) ( pp[ ) may now be 
expressed as 



IV. DISTRIBUTION OF THE PARTITION 
FUNCTION OF AN INFINITELY DEEP VALLEY 



Poo(p,u) = 



+ OC 



dye-wp^u) 



(30) 



dqR 00 (q)S ao (p,u,q) = fie- 13 " 



dss 



I, 



v ( Se ( , *s K (s) , 



I [s)Iu(s) 



Io(s) 



(31) 



In particular, the distribution of the energy u alone 
reads (taking into account the two sides y > and y < 0) 



Poo(u) = 2P oo (p = 0,w) = 2/3e 



- o/^-/ 31 ' / dss 





In \se ^ 



TO 



K o(se 2 )--rTTlo{se * 



(32) 



whereas the distribution of the position y alone has for 
Laplace transform 



Pooip) 



+ OG 



duP^p, u) = 2 



ds 



o sI Q {s)I v {s) 



dzzl v {z) (k (z)-^-Io(z) 



(33) 



The Laplace parameter p is present only through the in- 
dex v — of the Bessel function I v . 

The normalization to Pca(p = 0) _ 1 /2 for th e ha lf 
space can be checked using (A7) and ( AlC ). Using ( AH ), 
we may expand in v as follows 



Poo{p)= 

z -10 



dzz I $!A + o{S) 



h{z) 
KHz) 



1 4 /p /"» 
* / — / dzz 

2 3(3\JaJ I (z) 



(34) 

+ 0(p) (35) 



This shows that the probability distribution Poo{y) ex- 
hibits the power-law decay 



Poo(y), ,oc 



1 



2T 



K%{z) 
Io(z) 



(36) 



making all the integer moments diverging in the limit 
t — ► oo. The exponent 3/2 is of course related to the 



The partition function of the valley can be decomposed 
as the sum over two independent half-valleys 



dxe-WiW + / dxe-W^ (37) 
Jo 



Using the result ( |2q ) for Rr(q) (|2l|), we obtain that its 
probability distribution has for Laplace transform 



n+OO -, 

/ dZe-^Voo (Z) = Rl(q) = 

J 



F) 



(38) 



It is convenient to introduce the rescaled partition func- 
tion 



where 



Zoo 

AT 2 



(39) 



(40) 



represents the thermal length associated to the typi- 
cal scale of the extension of the Boltzmann distribution 
e -pu(x) j n a B r0WIuan well U(x) ~ a^fx. 

The probability distribution p{z) of the dimensionless 
partition function z has for Laplace transform 



dze- qz p(z) = - 



1 



o 2 (V?) 



(41) 



The series expansion (A3) shows that all the positive mo- 



ments are finite. The leading behavior at large z is indeed 
given via Laplace inversion by the first pole q\ = —s\ on 
the negative real axis, where s% is the first zero of the 
Bessel function Jq(s\) = : 



ze 



p{z) 

z — 

The behavior of fl4l| ) at large q 



dze~ qz p(z) ~ 27r^/ge 

g — ► oo 



-2^9 



(42) 



(43) 



leads to the following essential singularity at small z 

z^0 Z°'^ 



G 



V. LOCALIZATION PARAMETERS 



To compute the localizations parameters Ifc(oo) (17) 
we proceed along the same lines 



ife(oo) 



-Hoc. 



dy 



m Jo 



+00 



dqq^R^q) 



lim ( e ■"f" i v/e 
r^oo 



(45) 



r(*0 7 



+00 



fc-i 



i2oo(g) lim 7V"(r) 



r r 1 - 

j due- k ^ u lim ^F [0!r] (T - e, %)F [0;r] («, 0| 



e) (46) 



Using again the result (BIS) for -Ffo.r] we finally get 
fe-i 



Y k (oo)-- 



r(fc) v 4 
V^r 2 (fc) //3 2 



2fc-l 



2r(fc + i) V 4 



{) 

fc-1 



r 3 (fc) 
r(2fc) 



CSV)* -1 (47) 



where we have used again (|A7|) and ( |AT0| ). 

The increase of 1^. at large k is a consequence of the 
average over the disorder, and can be understood by con- 
sidering the averaged probability 1& to have k particles 
at the minimum of the valley (instead of at the same 
place but anywhere in the valley for Y k ) , which is exactly 
given by the negative moment of order k of the partition 
function Z x (37) 



2 



+00 



r(fc) \i T 



dZP^Z)^ 



+00 



ds- 



ISM 



(48) 



For large k, the dominant behavior comes from the the 
small z behavior ( ^3J ) of the probability distribution, 
which yields 



ft 



k — too 



r(k + -) 



(49) 



This shows that the behavior at large k of the average Y k 
is dominated by very steep valleys having a small parti- 
tion function z. 

Here we need to make some comments about the rela- 
tion with the discrete Sinai model with lattice constant a. 



For most quantities, the results obtained for the contin- 
uum version correspond to the universal limit where the 
lattice constant a is very small as compared to the ther- 
mal length It ~ l/(a/3 2 ) (40) representing the typical 
extension of the Boltzmann distribution in a Brownian 
well. For instance, this is the case for the probability 
distribution Poo(y,u) of the thermal packet and for the 
correlation function C(/,oo). For the localizations pa- 
rameters Yfe however, our result indicates that the di- 
mensionless y k parameters of discrete models will behave 



as 



yf screte (oo) 



V^r 2 (fc) 
2r(fc + i) \i T 



fc-i 



(50) 



when k is fixed, in the limit where (a /It) is small. But for 
fixed {a /It), there is a maximal value k max beyond which 
the result above doesn't apply anymore and is replaced 
by a non-universal behavior. Indeed, the discrete y k are 
by definition smaller than 1. And as explained above, 
the large k behavior of ( |47| ) is related to very steep val- 
leys having a small partition function z, i.e. involving 
a small number of sites in discrete models, so that all 
details of the model will be important to determine the 
large k behavior of yf lscrete (oo). 

Note also that in the opposite regime where the lattice 
spacing a is not negligeable with the thermal length It 
(i.e. a ~ It or a > It), there is only a small number 
of sites that are really important around the minimum 
of the valley so that the discreteness and details of the 
model will again be very important. For instance, the be- 
haviors of the localization parameters Yfc(oo) at zero tem- 
perature are highly non-universal and depend on many 
details : in (9) , the binary distribution of the random 
forces induces a lot of minima degeneracies separated by 
barrier of two bonds which can always be passed even in 
the the limit of zero-temperature (because the particle 
is not allowed to remain on the same site between t and 
t + 1). Assuming that all degenerate minima have the 
same weight, the value of 12(00) at T = is found to 
be (ln2)/2 M, instead of the 1 one would expect if there 
were no residual fluctuations at T — around a single 
minimum. 



VI. CORRELATION FUNCTION 

To compute the correlation function C{1, 00) ([l^), we 
decompose it into 



C(/,oo)= lim 2 / dqq / dy ( e'^^'e 



JO 



r^oo 



lim 4 / dqq ( e 



{U2} 



dy ( e 



, / e -0U 2 (l-y) e ~ q ff dxe-^2M \ 

{Ui} \ I {U2} 

iW \ 

-0(7x(2/)-/3i7 1 (3/+i) e -g/ o ' r dxe~^ u ^ \ 



{Ui} 



(51) 

(52) 
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which yields in Laplace transform with respect to I 



C{p, oo ) 



/•+0O rOO / 

/ dle~ pl C(l, oo) = 2 / <£qg lim / 
Jo Jo r ^°° \Jo 



r ,.+00 

J dye~ py S r (y,u,q) 



1 



(53) 



+4/ dqqR^q) lim Af(T) I du^-^j du 2 e-^F [0 , v] (u 2 ,p\ui) Hm -^F [0 , r] (T - e, 0|u 2 )F [0jr] (m, 0|e) (54) 



Using the previous result fl27|) and the expression (BIS 
for J^[o,r]j we finally get after simplifications 



r+00 

C(p,oa) = 8 dz 1 z 1 I ly (z 1 )K {z 1 ) 
Jo 

+00 

dz 2 z 2 K u (z 2 )K (z 2 ) (55) 



where again p only appears in the index v (|29| ) of Bessel 
functions. 



Expansion in p yields (A16) 



c(p,oc) = i-|y^+o(p) 



(56) 



This small-p behavior shows that C(Z, 00) presents at 
large distance the same power-law decay with exponent 
(3/2) as the probability distribution P<x>{y) (|36"|) 



(57) 



making again all the integer moments infinite. As ex- 
plained after equation (|3^), this long-range algebraic de- 
cay of the mean correlation function comes from rare con- 
figurations of the disorder, and is very different from the 
decay as e~' 3v ^ characterizing the typical correlations. 
This is thus an explicit example of the important differ- 
ences that exist in disordered systems between typical 



and mean correlations 1 1 1 



VII. COMPARISON WITH EQUILIBRIUM 
FUNCTIONS IN LARGE SYSTEMS 



In this section, we consider the Boltzmann equilibrium 
( |lO| ) on a finite system of length L to see if, in the ther- 
modynamic limit L — + 00, we recover the same properties 
for the thermal packet as in Sinai diffusion on the infi- 
nite line. Some statistical properties of the Boltzmann 
distribution ([h]) have already been studied in 
where in particular the decay of correlations was shown 
to be algebraic with the exponent 3/2, for the same rea- 
son as discussed above after the equation (|3^). The Y^ q 
have already been computed in |D| for free and periodic 
boundary conditions, but they were found to be very dif- 
ferent even in the thermodynamic limit L — > 00, whereas 
it is usually expected that for physical quantities that 



remain finite in this limit, differences should vanish. In 
the following, we compute the equilibrium functions for 
both boundary conditions and find that they coincide 
with each other, (i.e. the result (17) of (jl3| is erroneous). 
The thermodynamic limit is thus well-defined and inde- 
pendent of the boundary conditions. 



A. Y^ q parameters at equilibrium 



The Yi, q parameters for free boundary conditions may 



be rewritten in terms of path-integrals as 

rL 

rfree/ T \ _ 



1 i ■ 



dx- 



e -kpU(x) 







e -PU{y) 



T(k) 



dqq K ~ L / dxe- k P u We-<>So<lye- 0U(v) 
Jo 

+OO r-L /"+00 /> + 00 

dqq k ~ 1 I dx I du 



r(fc) Jo jo 

e~ kl3u G(u L ,L~- x\u)G(u,x\0) 
where the basic path-integral G is 



du 



(58) 
(59) 

(60) 



G(u,l\u ) 

U(l)=u 



U(0)=u 



VU(y)e-&tiM%)-9l? 



-pu(<e) 



(61) 



It is analogous to the path-integral (|23|) except that 
here there are no boundary conditions at U = and 
U = r, and the variable u is in ] — 00, +oo[. From a 
technical point of view, we mention here that contrary 
to the previous works ^-^] which expand the path- 
integral ( pl| ) upon eigenstates of the the Hamiltonian 
H = — + qe~^ u , here we have chosen to work in 
Laplace transform with resp ect to th e le ngth I to have 
a more compact result (C2). Using (C2), we obtain in 
Laplace with respect to L 



-Crfree/ \ 



+00 



dL, 



1 u 







r(fc) 



4 
ds 



k-2 



+oc 



A) / <> 1 - 



ds 



KM 



(62) 
(63) 
(64) 



with 
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2 [uj 



The thermodynamic limit L — » oo is obtained as 



lim Y^ ree (L) = limM^^M) 

L — ► oo (jj — *U 



rfree. 



2 / £ 

r(fc) ^ 4 



2r(fc + i) V 4 



dzz^^K^z) 



fc-i 



(65) 

(66) 
(67) 
(68) 



in agreement with the equation (19) of fl3[| . 

We now consider periodic boundary conditions, and 
indicate the modifications that appear. Taking into 
account that the probability to have U(L) — £7(0) is 
(l/v / 47TCTi), we have in terms of the path-integral ( |6l| ) 



-^periodic ^J^^ 



due 



-kpu 



-m-Jo dqq L dx 

G{0,L-x\u)G{u,x\0) 



(69) 



So here, it is simpler to compute the following Laplace 
transform 



+00 



dLe~ uL 


' xrperiodic 
k 


Vl 




r+oo 


/ dzz 
Jo 



(L) 



(70) 



2fc-l 



80F 

JO s J z s 

The thermodynamic limit L — > oo is then obtained as 

(72) 

v^r 2 (fc) //5 2 



lira F pe " 0<lk 'n - 

L — >oo 



(L) = hm — =2/ fc (w) 
fe-i 



2r(fc 



(73) 



contrary to the erroneous result in equation (17) of |[13| . 
Our result thus coincides with (|6^) concerning the equi- 
librium with free boundary conditions and with ( |4^ ) con- 
cerning the localizations parameters for Sinai diffusion on 
the infinite-line. 



B. Two-point correlation C eq (l) at equilibrium 

Similarly, the two-point correlation for free boundary 
conditions may be expressed in terms of the path-integral 



pL — l /*+oo pL — l f-\-oo 

C f L ree {l) = 2 dx[p e L q {x)p e L q (x + l)}=2 / dqq dx dqqe-PUW-PUV+^e-itfdye-' 317 ^ 

Jo " " Jo Jo Jo 

pL—l r+oo r+oo r+oo r+oo 

= 2 1 dx dqq dux du 2 e' l3ui ~ l3u2 du L G(u L , L - x - l\u 2 )G{u 2 , l\ui)G(ui,x\0) 

J J J — oo J — oo J —oo 



so that in double Laplace with respect to I and L we get 

c -foo pL 



p-\-oo rL oo /*+oo p 

Cf ree (p,Lu)= dLe~" L dle-P l C f L ree (l) = — dz 1 z 1 I v ,(z 1 )[K f ,(z 1 ) 

Jo Jo P ° Jo Jo 



ds 

(zi)[K^ Zl ) I -Us) 



+IJz x ) / —K^s)} / dz 2 z 2 K v ,[z 2 ) 



Z2 ds f°° ds 

KM I -I,(s)+I,(z 2 ) -K,{s) 



(74) 
(75) 

(76) 
(77) 



P 

is obtained as 

V/ree 



with v' = % \ . The thermodynamic limit L — ► oo 



lim C 

L — ► oo 



L (?) 



lim 

L — >oo 



dle- pl C L (l) 



(78) 



+00 



\im(ujC free (p,Lj)) = 8 / dziZiI v {zi)Ka{zi) (79) 



dz 2 Z 2 K v {z 2 )K Q {z 2 ) 



(80) 



pendix (^|) that periodic boundary conditions also yield 
the same result (jS0|). 



C. Probability distribution of the partition function 

The probability distribution of the partition function 

r-L 



Z, 



and thus coincides with the result ( poj ) for correlation of 
two particles of the thermal packet in Sinai diff usion on 
the infinite line. It is shown in equation (Cll) of Ap- 



(81) 



has already been computed in p4j , but to compare with 
our result (5^) for the infinitely deep valley, we need to 
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consider the modified partition function 



(82) 







where U m i n is the minimum of U(x) for < x < L. Us- 
ing the notations of (|23|), the Laplace transform of the 
probability distribution Pl(Z) can be expressed in terms 
of path-integrals as 



P L ( 9 )= / dZe-* z P L (Z) = du du L / dx 

JO J-OO Jtln JO 



Uo 

e Ju(o)=o 



U(L)=u L 



E>Z7(x)e w Jo <*«l*J 6 [uo , +oo] {C7(a;)} 



'!7(xo)=«o+e 

The change of functional ^(ar) = U(x) — uq leads to 

/•+oo /"-t-oo />L -i pV(xq)— € 2 

A,(?)= / d«o / d«i / dx Q \im- VU(x)e-^^ L ^m-^ L d X e-^^ Q {u{x)} 



(83) 



JO Jo 

V(i)=«£ 



e^O 



V(0)=«o 



V(x )=e 



T>V(x)e 4 " J " V^; «Jo e[ ,+oo]{C^(^)} 



= (7 / rfu / d«L / rf^o hm — F[o !+oo] (e,a; |w )i 7 '[o, + oo](wL,i- a^o|e) 
Jo Jo Jo e ~*° e 



(84) 
(85) 



whic h yie lds in Laplace with respect to the length L, 
using (gig) 



+ 00 



dLe-" L P L (q) 



1 /"+ 00 
lim - / duF [0 +oo] (e, 



7 C / — 



M /3V a 



(86) 
(87) 

(88) 



where has been defined in (65). Again, the thermody- 
namic limit L — > oo is obtained as 



D. Conclusion 

The conclusion of this section is that the statistical 
properties of the thermal packet for the Sinai diffusion 
in the infinite time limit exactly coincide with the ther- 
modynamic limit L — > oo of the statistical properties 
concerning independent particles at Boltzmann equilib- 
rium in a sample of size L, with no dependence on the 
boundary conditions. 



VIII. LOCALIZATION PROPERTIES OF 
FOKKER-PLANCK EIGENFUNCTIONS 



Poo {q)= lim ( u I dLe-" L P L (q) 



T 2 2 fq 
J UVa 



(89) 
(90) 



which coincides with the result ( |38| ) for the infinitely- 
deep valley. It is of course straightforward to gener- 
alize this computation to obtain that the joint distri- 
bution P QO {y 1 u) ( |3ll) also coincides with the thermody- 
namic limit L — ► oo of the finite-size joint distribution of 



j) and (u — U m in) where U m i n is the minimum 



of U (x) for < x < L and where x r 
this minimum. 



is the position of 



As discussed in Section ||, the approximation (15) be- 
comes asymptotically exact in the limit r — > oo. It is 
thus interesting to explore the consequences of this ap- 
proximation for the eigenfunctions of the Fokker-Planck 
operator. 



A. Recall of some exact results 

To discuss the properties of eigenvalues and eigenfunc- 
tions of the Fokker-Planck equation (Q), it is more con- 
venient to use the well-known transformation into an 
imaginary-time Schrodinger equation via the introduc- 
tion of the Green function 

G{x, t\x ,0) = e^ u ^~ u ^P{x, t\x , 0) (91) 

which satisfies 
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d t G(x, t\x , 0) = ~H s G{x, t\x , 0) 



H s = -Tdl + { — U'(xf 



1 



U"{x)) 



(92) 
(93) 



with the initial condition G(x,t\xo,0) — > 8{x — xq)- 
This is the standard form for the Schrodinger operator 
Hs associated to a diffusion process. It can be factor- 
ized as H s = TQ^Q with Q = d x + U'{x)/{2T) and 
= —d x + U'(x)/{2T), and has thus a real positive 
spectrum. We consider the case of a very large but finite 
system where the spectrum of energies E n is discrete. 
The Fokker Planck operator Hpp is non hermitian but 
has the same real positive spectrum, with right and left 
eigenfunctions ^„(x) an d ®n( x ) associated to E n . They 
are related to the eigenfunctions ip n {x) of the Schrodinger 
operator by 



,U(x)/(2T) 



(94) 
(95) 



The expansion upon Fokker-Planck eigenfunctions now 
reads 



P(xt\x 0) = y$«(x)^(x )e 



„-E n t 



(96) 



The ground-state n ~ has of course for energy Eq = 
and corresponds to the relaxation towards Boltzmann 
equilibrium so that the left and right ground state eigen- 
vectors are simply given by : 



$o (*)= 




(97) 
(98) 



where Z to t = J dxe U ^ X ^ T is the normalization over the 
finite large system. 



B. Construction of an orthonormalized set of 
eigenfunctions for the approximate dynamics 

In this section, we consider the approximation (|l^) for 
the dynamics as defining a new dynamics denoted by tilde 

P(xt\x o 0) = V J-e-WWevt (x)6 Vr (x ) (99) 
and try to expand it upon a basis of eigenfunctions as in 



P(xf|x 0) = ^$«(x)^(x )e-^ 



(100) 



At time t, the states in ( Jl 0C| ) with energies E n > 1/t 
are negligeable in the sum, and correspond in the RSRG 
picture to bonds that have been already decim ated . The 
state n will become negligeable in the sum (10C) at a 
time t n ~ l/E n , and this disparition from the sum cor- 
responds in the RSRG picture to a decimation at scale 



r„ = T\nt n — —TlnE n . The low-lying energies E n are 
thus exactly determined by the large RG scales T n at 
which decimations occur in the system. Of course in the 
real system, in the case of near degeneracies of neighbor- 
ing bonds, slight shifts in these levels will occur. 

To determine now the eigenfunctions, we consider what 
happens upon this decimation at T n . We may consider 
that the time exponential factor associated with level n 



has changed from e 



1 to e 



-E n t+ 



while none 



of the others decaying exponentials in (100) have changed 
(since time scales are well separated). The difference 
P(xt~ \x 0) - P(xt+\x o 0) is thus equal to ^(x)^(x ) 
On the other hand from the RSRG point of view, what 
happens is simply a decimation where two valleys V\ and 
V2 have merged into a single one V . Thus we get using 



6 Vl {x)6 Vl (x ) e V2 (x)9 V2 (x ) 



(9 Vl (x) 



1v 2 (x))((0vi(x o ) + Ov^Xq)) 



z 



Vi 



(101) 
(102) 

(103) 



One sees that indeed the right hand side has the 
nice property to factorize into a function of x and 
one of xq and thus one can determine unambiguously 
&n(x) and $^(xo) by fixing the constant using $^(x) = 
e _!7 ( x )$^(x) from (|95|), which leads to the eigenset for 
n > 1: 



*n(*) = \ f V +y (^0 V1 (X)-^9 V2 (X)) (104) 



(105) 



One can check on ([105|) that the eigenfunctions have 
all the correct normalization and orthogonality proper- 
ties. First one has for n > 1 



dx$%(x) = 



(106) 



This insures the normalization of the probability distri- 
bution P(xt\xoO) for all t and xo as it should. 



dxP(xt\x o 0) = I dx^(x)^(x ) = 1 (107) 



Second one finds the correct normalization 
J dx^{x)^(x) = J dxe- u W[$%(x)] 7 = 1 (108) 

Furthermore, one can also check that the set of wave- 
functions exactly forms an orthonormalized set : 



(109) 
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So the disorder-dependent form (105) for the eigen- 
functions have all t he good properties to represent via 
the expansion (100) the dynamics defined in (|99"|). One 
may even define the following effective Fokker-Planck op- 
erator 



H FP = J2E n \^><^\ 

n 

as an approximation to Hpp (§)■ 



(110) 



C. Qualitative properties of eigenfunctions 

We can now discuss the typical shape of an eigenstate 
and compare with the qualitative features discussed in 
the Appendix B of Q] for the related model of one- 
dimensional random-hopping hamiltonian for fermions 
(except that here there is no particle-hole symmetry and 
there is no need to disti nguish even and odd sites) : 

(i) an eigenstate (105) has two peaks corresponding to 
the minima of valleys V\ and V^- The eigenstate has a 
significant value in finite regions of order It ~ ( ffO| ) 
around these two peaks. 

(ii) these two peaks are separated by a distance of order 

r 2 ~(in£:) 2 . 

(iii) Away from one of these bumps but within the 
valley, ie. for distance r < T 2 , the decay of the wave- 
function is governed by the Boltzmann factor 



-/3U(r) 



which typically behaves as e" 



Cv ^. In particu- 
T 2 , this gives 



lar, on the edges of the valleys where r 
an amplitude of order e _c r . 

(iv) Beyond t he in volved two valleys, within our simple 
approximation (105) with theta- functions on the edges of 
valleys, the eigenstate is simply zero. So to estimate the 
decay of the eigenstate for distances r > T 2 , we must take 
into account as in the Appendix B of Jl4| that the two 
points are typically separated by a number ^ of valleys, 
and that the overlap between two neighboring valleys is 
not exactly zero but of order e _c r . And thus a pertur- 
bation theory yields that the decay for distances r > T 2 
behaves as e~ c r. We thus recover that the localiza- 
tion length in the sense of asymptotic exponential decay 
of the associated quantum wavefunction thus behaves as 
X(E) ~ r ~ (— TlnE) whereas the typical extension of 
an eigenstate behaves as the typical distance between the 
two peaks which behaves as 1(E) ~ T 2 ~ (— TlnE) 2 and 
which directly gives the low-energy behavior of the inte- 
grated density of states N(E) ~ l/l(E) ~ l/(-Tln£') 2 
||. We note that the length 1(E) is the one which appear 
in the averaged Green function ( |9l| ) as computed in ]22] ] 
for the lattice fermion model with random hopping, and 
computed in Q| for the present prob lem. 

So the disorder-dependent form (105) for the eigen- 
functions, even if only approximate near the edges of the 
valleys, give a good insight into the properties of low- 
energy eigenfunctions. For a study of systematic correc- 
tions to this approximation, and further results on the 



statistics of wavefunctions far in the tails, we refer the 
reader to D.S. Fisher §1. 



IX. PROPERTIES OF METASTABLE STATES IN 
SINAI DIFFUSION 

In the slowly relaxing systems such as glasses, gran- 
ular media or disordered spin models, the notion of 
metastable states is an important and ubiquitous con- 
cept. In particular, the Edwards ergodicity conjecture 
|17| consisting in computing dynamic quantities by tak- 
ing flat averages over metastable states has given rise to 
a lot of recent studies fl^ , |IgH 2l"|| . However, as discussed 
in details in p"5| , truly metastable stables only exist in 
mean-field approximations or in the zero-temperature 
limit, and thus to use this concept at finite tempera- 
ture in finite-dimensional systems, one needs to consider 
metastable states with finite lifetimes. In this context we 
find instructive to consider what happens in Sinai diffu- 
sion from the point of view of metastable states and to 
compare with the general theory of glassy systems. 



A. Identification of the set of metastable states at a 
given time t 

In the general theory of glassy systems, it is natu- 
ral to separate the dynamics in two parts : there are 
'fast' degrees of freedom which rapidly reach local quasi- 
equilibrium plus a 'slow' non-equilibrium part. At a given 
long-time t, the fast motion covers a region of phase-space 
which can be defined as a metastable state associated to 
time t ]l5| ]. In Sinai diffusion, this picture directly ap- 
plies. The formula ( |l5|) leads to a very direct identifi- 
cation of metastable states : at time t, the metastable 
states (i) are given by all the renormalized valleys V-ff* 
that exist at scale T = Tint. Moreover, the formula 
( |l5| ) exactly corresponds to the construction of Gaveau, 
Schulman and Lesne fl6|| as summarized in [p"5"|| , where 
the evolution operator e~ tHpp is replaced by a projec- 
tor onto the subspace of states (i) having an eigenvalue 

E n < 1/t 



- tH FP ~ ^ I Pj >< 



(111) 



with the following identifications : the right eigenvectors 
read 

p i( x ) = - — i—mw) e{x G F r w ) (112) 



and thus they have indeed the properties of being posi- 
tive, normalized and not-zero in non-overlapping regions 
of space. The left eigenvectors read 



(113) 
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They are thus indeed one within the support of Pi (x) , and 
zero everywhere else. Moreover, the fast degrees of free- 
dom have indeed converged towards a local Boltzmann 
equilibrium within the meta stable state represented by 
the renormalized valley (112). 



B. One-time dynamical quantities as averages over 
the set of metastable states at a given time t 

We now consider the question of Edwards conjecture. 
In Sinai diffusion, within the RSRG approach (?]] , all one- 
time quantities are effectively computed by averages over 
all renormalized valleys, with a measure which is not a 
priori flat, but which depends on the computed quantity. 
In particular, for a uniform initial condition, the spatial 
length of a renormalized valley exactly represents the size 
of the basin of attraction of this renormalized valley. So 
the rescaled size A = j& = (Tint) 1 °^ * ne basins of at- 
traction of the metastable states is a random variable 
distributed with the probability distribution 



C. Hierarchical organization of metastable states 
associated to different times 

In the case of Sinai diffusion, one knows much more 
than the one-time description in terms of metastable 
states. Indeed, the expressions ( |105| ) for the Fokker- 
Planck eigenvectors contain in fact all the information on 
the changes in time of the set of metastable states. They 
have an obvious hierarchical organization : two eigen- 
states can be either disjoint or nested (ie. one is included 
in another) . And the evolution in time is described by the 
RSRG procedure. This is why two-time aging quantities 
can be computed within the RSRG method f7j,p6|. 



D(\) = P*(.)* X P*(.) 



(114) 



where P*(X) is the distribution of the rescaled length of 
bonds M. More explicitly, its Laplace transform reads 



D(p)= P*{p) 



1 



cosh y/p 



(115) 



and inversion yields |2[ 



n— — og 



2tt z [ n + - | A - 1 



+00 



7^ £ ( - ir+1 

v m— — oo 



me 



e -^ 2 (n+i) 2 (n6 ) 
(117) 



Since Edwards conjecture is usually based on the as- 
sumption that all the basins of attraction of the various 
metastable states have the same size |l5| , it is of course a 
very strong hypothesis that cannot be true in general but 
only for special systems with special dynamics [ p"8pl| ] . In 
Sinai diffusion, the metastable states have not the same 
size and thus Edwards conjecture cannot be satisfied in 
general. 

However, for special quantities, Edwards conjecture 
can be recovered. For instance, in this paper, we have 
computed the probability distribution Pco(jJ,u) of the 
thermal packet, the localization parameters Yfc(oo) (j|) 
and the correlation function C(Z,oo) (|]) as flat averages 
over infinitely deep wells. This is because in the limit 
r — > oo, only the lower part of the Brownian valley is 
important, and thus there is no dependence on the size 
of the valley for these quantities concerning the thermal 
packet. 
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APPENDIX A: USEFUL PROPERTIES OF 
BESSEL FUNCTIONS 



APPENDIX B: EXPLICIT COMPUTATIONS 
FOR INFINITELY-DEEP VALLEYS 



The Bessel functions I v {z) and K v {z) are two linearly 
independent solutions of 



z 2 f"(z) + zf'(z)-(z 2 + v 2 )f(z) = 
Series expansion in z reads 



+00 



1 {Z\ 2k 



+00 



k=0 1 '' 



^ ip(k + 1) /z\ 2k 

Al- 



and 



+00 



K v {z 



1 



k\T{v + k + 1) V2 

TV I- V {z) - I v {z) 

2 sin 7rzy 



The wronskian property 



I v {z)K' v {z)-I' v {z)K v {z) = — 

z 



gives 



+00 



ds- 



+00 



ds- 



+00 



ds- 



1 


'Ms)' 


"f 








h(s) _ 




Hz) 


K (s) 1 


\K%{ S 


)" 


+ °° * 2 (*) 


s/ 3 ( S ) 2 


U 2 (*)_ 


, ~ 2ig(*) 


Ki(s) 1 




)" 


+ °° * 3 (*) 


S / 4 ( S ) " 3 


Uo 3 (*)_ 


3/ 3 (z) 



Another useful integral is for any k > 



/o 



4r(fc + i; 



The equation (Al) leads to 

•2 



d_ 

ds 
<]_ 
ds 
d 
ds 



S ^(K 2 (s)-[K' (s)f 



= sK 2 (s) 



-(I (s)K ( S )-I' (s)K' (s)) 



= sI (s)K {s) 

Differentiation with respect to order 

I v {z)= I {z) - vK (z) + 0{v 2 ) 
K u {z)=K Q {z) + 0{v 2 ) 



(Al) 

(A2) 
(A3) 

(A4) 
(A5) 

(A6) 

(A7) 
(A8) 
(A9) 

(A10) 

(AH) 
(A12) 

(A13) 
(A14) 



(A15) 
(A16) 



1. Explicit expression for K r](u, x\uo) 



The path-integral FmriOw, x\uq) defined in ( |23| ) can be 
obtained via the Feynman-Kac formula. Indeed, as a 
function of the variables (u, x), it satisfies the imaginary- 
time Schrodinger equation 



d x F = ad 2 u F - qe~ f3u F 
with the initial condition at x = 
F(u,0) =S{u-u ) 



(Bl) 



(B2) 



and the absorbing boundary conditions at u = and 

u = r 



F(0,x) = 
F(T,x) = 

The Laplace transform with respect to x 
r+00 

F(u,p) = / dxe- px F(u, x) 
Jo 

satisfies the system 

ad 2 u F - qe- f3u F - pF = -8{u - u ) 
F(0,p) = 
F(T,p) = 

Let us introduce two linearly independent solutions 
<pi(u,p) and <p2(u,p) of the equation: 



(B3) 
(B4) 



(B5) 



(B6) 



ad 2 u F - qe- pu F - pF = 
In terms of the new variables 

2 fq 
p V cr 
= 2 fp 
0\ a 

the equation becomes 

z 2 F"{z) + zF'{z) - {z 2 + v 2 )F{z) = 



(B7) 

(B8) 
(B9) 

(BIO) 



Two linearly independent solutions of this equation are 
the Bessel functions I v {z) and K v {z) so we choose 



k(u,p) = I v {z) = I 2 fT\\\\-e- 
?v <r \p V a 



<fc(u,p) = K v {z) = K 2 



0y a 



(Bll) 
(B12) 
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and the wronskian of these two solutions reads 

to = 0i(u)^(u) = § (B13) 
We now introduce the function 
E(u,v,p) = 

— {<l>i(u,p)fo(v,p) - <p 2 (u,p)(j>i(v,p)) (B14) 
w 



and express two solutions of the equation (B7) which 
vanish respectively at and T as: 

$_(u,p) = J5(0,u,p) (B15) 
<Z> + (u,p)=E(u,T,p) (B16) 

Their wronskian reads : 

W(p) = *'_(«, p)*+(u,p) - $_(w,p)$' + («,P) 

= E(Q,r,p) (B17) 



The solution of the system (B6) now reads : 



F[o,r]( u >P\ u o) = 

(min(u, M ),p)$ + (max(u, u ),p) _ 
aW(p) 

E(0, min(it, uo),p)E(max.(u, uq), T,p) 

<r£(o,r )P ) 



2. Explicit expression for Roo(q) 



(B18) 



Specializing ( B18| ) to the case uq = e and it = T — e 
and expanding 



E(0,e,p)=e + O(e 2 ) 
E(T-e,T,p)=e + 0(e 2 ) 



we get 



!^^ F [°.r]( r - £ 'Pl £ ) 



<rS(0,r,p) 



Then @ yields 
Rr(q)= Af(T) lim -^F [0 ,r] (T - e,p = 0|e) 

Af(r) 



(B19) 
(B20) 



(B21) 

(B22) 
(B23) 



aE(0,T,p = 0) 

Using the expansions ( [A3] ) of Bessel functions, the nor- 
malization condition Rr(q — * 0) = 1 determines the nor- 
malization 



W(r) = ar 

so that the final result reads 



Rr(g) 



2 



Jo(s)^ (se-^) - K (s)I (se-^) 



(B24) 
(B25) 



with s = | J?-. 



3. Explicit expression for Sr(y,u,q) 

We now turn to the evaluation of ( p6| ) in Laplace with 
respect to y 



Sr{p,u,q)= / dye pv S T (y,u,q) 



(B26) 



N{T)e~^ lim ( 4riWr - e,0|u)F [o ,n(«,p|e) ) (B27) 



Using 



r f 1 * ( \\\ %M 



and 



jg& {- e F ^ r -^) = aE (o,r,Q) 



(B28) 



(B29) 



we get 
Sr(p,u,q)-- 



_I_ e -,. E( „,„, 0) |^ (mo) 



In the limit T — > oo, using the expansions (A5) we get 



lim 



E(u,T,p)\ '»\JS 



and 



r^oo vs(o,r,p) 



lim 



e 2 



(B31) 



(B32) 



so that we obtain the formula (^]) given in the text 



APPENDIX C: COMPUTATION OF 
EQUILIBRIUM FUNCTIONS IN LARGE 
SYSTEMS 



1. Basic path-integral 

The Laplace transform with respect to I of the basic 
path-integral ( fill ) satisfies the same equation as (B6) but 
the boundary condit ions are now at u — ► ± oo. S o using 
again the notations (BE) the solution reads ( B18 ) 



G{u,p\u ) 



j3a 



(CI) 



(u,u )\ ^ QyZ e -#max(«,» )^ ( C2 ) 
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2. Two-point correlation C cq (l) for periodic 
boundary conditions 

For periodic boundary conditions, the two-point cor- 
relation may be expressed in terms of the path-integral 
©as 



Cf H °* c (0 = 2Vi^I / dx dqq 
Jo Jo 

"+00 /*+oo 



dme-^ 1 / du 2 e~ 13112 

J — OO 

G(0, L — x — l\u 2 )G(u 2 , l\ui)G(ui,x\Q) 
which yields in Laplace transform 



dLe~ uL / dle 
^0 



c(p,w) 



po- Jo s J 
[K^I^eiz! -s) + In(zi)Kn(s)0(a - z x )\ 



+ OO 



d,Z2Z2K v t{z2) 



(C3) 

(C4) 
(C5) 



(C6) 

(C7) 
(C8) 
(C9) 



[K^(z 2 )I fl (s)6(z 2 -s) + I^z 2 )K^s)6{s - z 2 )] (CIO) 



with 11 = -s\f^- and v — 

/J V cr p 



1 1 "' - -y 2 ^. The thermodynamic 

limit L — » 00 is obtained as 



Cr*(p)=lim(^a(p )W )) 



+00 



8 / d2i2i/„(2i)i^ (^i) 



^2^2^(22)^0(22) 



(Cll) 



APPENDIX D: RELATION WITH GOLOSOV 
THEOREM 

As explained in the introduction, the theorem of 
Golosov Q states that -Poo(y) is equal to Pg{v) given 
by equation (||) . The purpose of this appendix is to show 
this formula of Golosov gives the same result as the equa- 
tion (|3^) obtained in the text. We first rewrite (|^) as 



L — >oo 



P G (y) = lim(/ dq Il e -«S L ^ {t) 



(Dl) 



»„> < D2 > 



Since the Bessel processes {r(t)} and {p{t)} may be 
considered as radial parts of free three-dimensional Brow- 
nian motion, the Feynman-Kac formula yields 



q Jf die-vow 



e *jo 
+00 



l-t-OO 

/ 4TrR 2 L dR L G p (R L ,L\0) (D3) 
Jo 

■Wile-^, 1 ^*')) ^ = j 4nR 2 L dR L 

/ AirR 2 dRG q (R L ,L-y\R)e- 0R G q {R,y\O) (D4) 
Jo 



where G q (R, Z|0) satisfies the imaginary-time Schrodinger 
equation 



diG q (R,l\Ro) = —H q G q (RJ\Ro) 



where 



a d 



H ^-R^dR{ K ~dR^ +qe 



R 2 — \ - :J! 



(D5) 



(D6) 



is the radial restriction of the corresponding 3D Hamil- 
tonian. We have moreover the initial condition 



G q (R,l^0\Ro) 



1 



AttR 2 



S(R-Ro) 



(D7) 



The identity (^im) = ^W), leads 

the change of function 



to 



g q {R, l\Ro) = 4<KRRoG q (R, l\R ) 



(D8) 



This new function g q satisfies the one-dimensional 
Schrodinger equation 



d l9q (R,l\Ro) = -h q9q {R,l\Ro) 
d 2 



h * = - a dR? +qe 



-pR 



(D9) 
(D10) 



on the semi-infinite line R > with the absorbing bound- 
ary condition at R = 



g q (R = 0,l\Ro) = 
and with the initial condition 

g q (RJ^O\R )^S(R-Ro) 



(Dll) 



(D12) 



By comparison with the equations (B1-B2-B4), we im- 
mediately obtain 

g q (R,l\R ) = hm F [0<r] {R, l\Ro) (D13) 

Since the function G q (R, Z 1 0) has to be obtained as the 
limit 



G q (R,l\0) = ^lim(- e g q (R,l\e] 
we have, using flBS|) and (|A5|) 



(D14) 



1G 



e H] o 



{p} 



\hn J dRRlimf^g q (R,L\e) 



lim 



(§-) 2 



I v (*/X)Jo V/3 V ^ 



(i V a 



(D15) 
(D16) 
(D17) 



in agreement with (|25|). 

Similarly, we compute the Laplace transform with re- 
spect to y 



1 



+00 



dye-™ ( e -M») e -« Jo + ~ *^' w ) (D18) 

/•OO />oo 

= lim[s/ dRi? / dR'gJR lS \R')e-P R ' (D19) 
Jo Jo 

hmQ§ ? (i?',f>|e))] (D20) 
which gives the same result as (p7[). 
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